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Introduction  and  objectives 

A  critical  requirement  for  the  development  of  low  cost,  wide-bandwidth 
telecommunications  equipment  is  the  close  integration  of  both  digital  and  analog  microelectronic 
components.  The  physical  interface  between  an  IC  and  its  environment  is  the  IC  package,  and  its 
performance  is  severely  tested  by  the  high  speed  and  high  frequencies  encountered  in  wide- 
bandwidth  systems.  Single  chip  mixeu  signal  ICs  that  combine  directly  both  high  frequency 
analog  and  high  speed  digital  sub-sections  will  require  proper  electromagnetic  understanding  of 
capacitive,  inductive,  and  radiative  coupling  between  components,  and  their  impact  on  high 
sensitivity  analog  sub-circuits. 

The  overall  objective  of  this  work  is  the  development  of  techniques  for  the  analysis  and 
design  of  mixed  signal  packages,  especially  the  impact  of  inductive  cross  talk  between  the  digital 
and  analog  sections  of  the  IC,  and  techniques  to  maintain  appropriate  RF  to  RF  signal  line 
isolation  in  low  cost  IC  packages.  In  wireless  personal  communication  services  (PCS)  units  there 
are  particularly  severe  constraints  on  package  options,  since  these  applications  are  typically  very 
cost  and  form-factor  sensitive.  Some  of  the  major  problems  induced  by  digital  to  analog  and 
analog  to  analog  cross  talk  are  listed  below: 

(1 )  high  speed  digital  clocks  cause  severe  interference  with  RF  or  IF  front  ends; 

(2)  in  digital  portables,  time-division-multiple-access  (TDMA)  may  be  used  and  power  on  / 
off  cycles  happen  fairly  frequently,  causing  additional  transient  noise  on  the  power  and 
ground  planes; 

(3)  in  frequency  division  duplex  systems,  high-power  transmit  signals  cause  interference  with 
weak  receive  signals  since  separation  by  filters  is  limited; 

(4)  the  leakage  of  the  amplifier  output  to  the  input  may  cause  the  amplifier  to  oscillate. 

Under  this  program,  we  are  developing  and  testing  electromagnetic  modeling  techniques 
that  can  capture  such  effects;  to  achieve  this  it  is  critical  that  the  electromagnetic  analysis  tools 
exclude  unnecessary  effects,  but  without  requiring  "expert"  intervention.  We  are  developing 
inductive  extraction  processes  that  use  ult  a-compact  equivalent  circuits  (consisting  of  frequency- 
independent  elements)  to  model  frequency  dependent  skin  and  proximity  effects.  Recent 
accomplishments  include  a  new  highly  efficient  simulation  method,  the  surface  ribbon  method 
(SRM),  that  allows  much  more  rapid  simulation  to  be  performed.  Typically,  computation  time 
can  be  reduced  by  about  a  factor  of  2000  compared  the  conventional  model.  We  are  currently 
preparing  user  software  for  transmission  line  impedance  modeling  based  on  the  SRM,  and  expect 
release  in  first  quarter  1997.  We  have  recently  derived  a  rigorous  definition  for  the  “effective 
internal  impedance”  used  with  the  surface  ribbon  method,  and  will  begin  investigating  its  use 
with  the  finite  element  method.  This  approach  should  allow  significant  improvement  in 
computation  time  for  FEM-based  electromagnetic  simulations  of  finite  conductivity  structures. 
We  are  also  studying  intemal-to-the-package  low  pass  filter  designs  for  analog  power  supply  to 
chip;  pad  /  pin  arrangements  for  reduced  digital  power  plane-to-analog  coupling,  for  RF 
controlled  impedance  connections,  and  for  increased  RF-to-RF  signal  line  isolation.  Prototypical 
wireless  communications  circuits  will  be  used  to  determine  the  impact  of  digital-to  analog 
interference,  and  validate  the  models  developed. 
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Accomplishments 

Effective  Internal  Impedance  Boundary  Condition  for  Finite  Conducttvty  Metals 

Impedance  boundary  conditions  (IBCs)  are  widely  used  in  scattering  problems,  eddy 
current  problems,  and  lossy  transmission  line  problems.  An  DBC  is  usually  adopted  to  simplify 
the  problems  by  eliminating  from  the  domain  to  be  solved,  for  example,  lossy  dielectrics,  multi¬ 
layered  coatings  on  conductors  in  scattering  problems,  or  the  lossy  conductors  in  eddy  current 
and  transmission  line  problems.  Solutions  can  then  be  obtained  by  applying  the  finite  element 
method  (FEM)  [1,  2],  the  boundary  element  method  (BEM)  [3,  4],  the  electric/magnetic  field 
integral  method  (MF1E/EFIE)  [5,  6],  the  finite  difference  time  domain  method  (FDTD)  [7,  8], 
etc..  In  each  case,  the  use  of  an  IBC  reduces  the  number  of  unknowns  and  can  substantially 
reduce  the  computation  time  required.  But,  in  general,  the  boundary  condition  must  be  known 
(at  least  approximately)  a  priori  on  the  surface,  and  the  accuracy  of  the  field  solution  is 
determined  by  the  accuracy  with  which  the  IBC  is  approximated. 

The  most  widely  used  IBC  is  the  standard  impedance  boundary  condition  (SIBC),  also 
called  the  Leontovich  boundary  condition  [9, 10],  that  was  originally  developed  for  use  when  the 
skin  depth  is  small  relative  to  other  dimensions  of  the  problem,  or  a  layer  is  thin  and  highly  lossy. 
In  general,  for  coated  dielectric  layers  on  scatterers,  the  reflection  of  the  electromagnetic  wave  at 
the  boundary  depends  on  the  angle  of  the  incident  wave,  and  hence  requires  higher  level  of  IBCs 
[10],  such  as  the  Tensor  Impedance  Boundary  Condition  (TTBC),  the  Higher  Order  Impedance 
Boundary  Conditions  (HOIBC),  etc.  For  eddy  current  and  transmission  line  problems,  the 
Leontovich  boundary  condition  has  been  used  to  model  lossy  conductors  at  high  frequency.  In 
this  paper,  an  efficient  approach  for  calculating  the  series  impedance  of  lossy  multi-conductor 
lines  from  DC  to  high  frequency  is  presented  that  use  the  effective  internal  impedance  (EH)  as  an 
impedance  boundary  condition.  The  advantage  of  the  EH  when  compared  to  the  SIBC  is  that 
even  for  strongly  coupled  transmission  lines  it  is  easily  approximated  from  low  frequency  (i.e., 
the  skin  depth  is  larger  than  the  cross-sectional  dimensions  of  the  conductors)  to  high  frequency 
(i.e.,  the  skin  depth  is  far  smaller  than  the  dimensions  of  the  conductors)  using  the  surface 
impedance  of  an  isolated  conductor.  As  an  example,  the  EH-based  method  and  the  SIBC 
combined  with  the  BEM  are  compared  in  the  case  of  twin  lead  conductors,  and  the 
appropriateness  of  the  EH  approach  is  shown  for  series  impedance  calculations  of  lossy  quasi- 
TEM  transmission  lines  from  DC  to  high  frequency. 

Schelkunoff  [11]  first  introduced  the  concept  of  surface  impedance  in  electromagnetics  in 
1934  for  the  analysis  of  coaxial  cables.  In  the  1940’s  Leontovich  [9]  as  well  as  many  other 
workers  studied  the  basic  properties  of  the  surface  impedance  for  a  semi-infinite  plane  of  an 
isotropic  linear  medium  and  on  a  conductor-backed  thin  lossy  dielectric  layer  where  a  plane  wave 
is  incident.  Senior  [12]  explained  in  detail  the  Leontovich  boundary  conditions  and  the  requisites 
to  be  satisfied. 

In  order  to  approximate  the  SIBC  for  multi-conductor  problems  the  relationship  between 
the  tangential  electric  and  magnetic  fields  at  any  point  on  the  boundary  must  be  (at  least 
approximately)  a  local  one,  depending  only  on  the  curvature  of  the  surface  and  the  local 
electromagnetic  properties  (i.e.,  |i,  e,  and  a)  of  the  bodies.  Hence  the  surface  impedance  must 
not  be  altered  by  the  global  geometry  of  the  conductor  or  the  existence  of  other  conductors  in  a 
problem.  This  constraint  will  hold  if  the  operating  frequency  is  “high”  (but  much  lower  than  the 
dielectric  relaxation  frequency),  i.e.,  when  the  curvature  radius  of  the  surface  is  larger  than  the 
skin  depth,  when  the  refractive  index  of  the  bodies  is  larger  than  that  of  the  external  medium,  and 
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when  the  dimensions  of  the  problem  are  smaller  than  the  wavelength.  In  such  circumstances,  the 
reflection  characteristic  is  independent  of  the  incident  angle.  For  geometries  having  curvature, 
Leontovich  introduced  a  first  order  curvature  correction  term  to  the  surface  impedance  for  small 
radii  of  curvature  and  Mitzner  [13]  later  refined  this. 

At  low  and  middle  frequency  where  the  skin  depth  is  comparable  to  or  larger  than  the 
dimensions  of  the  conductors,  the  surface  impedance  is  strictly  no  longer  a  local  property  and 
depends  on  the  global  geometries  of  the  conductors.  For  “thin-film”  transmission  lines  using 
metals  with  thicknesses  of  less  than  a  few  micrometers,  such  a  condition  may  hold  for 
frequencies  well  into  the  GHz  range.  Under  such  conditions  the  surface  impedance  when  other 
conductors  are  present  may  considerably  differ  from  the  surface  impedance  of  an  isolated 
conductor.  Hence,  it  is  difficult  to  know  or  approximate  a  priori  the  surface  impedance  in  lossy 
multi-conductor  systems  at  low  to  middle  frequencies.  This  limits  the  usefulness  of  the  SIBC  for 
calculating  the  series  impedance  of  lossy  quasi-TEM  transmission  lines  at  low  and  middle 
frequencies,  and  necessitates  complicated  models  for  the  IBC. 

It  would  be  useful  to  find  some  technique  in  which  surface  characteristics  can  be  more 
easily  approximated  than  the  standard  surface  impedance  boundary  condition.  Instead  of 
developing  complicated  SIBC  models  to  capture  the  coupling  and  non-localized  field  effects 
between  multiple  conductors  at  low  and  middle  frequency  it  is  possible  to  develop  a  formulation 
in  which  isolated  conductor  surface  impedance  is  a  useful  approximation.  In  this  approach,  the 
conductor  interior  is  not  excluded  from  the  domain  of  solution,  but  rather  the  conductor  interior 
is  replaced  by  the  exterior  medium  and  the  Green’s  function  of  the  exterior  medium  is  used  also 
for  the  interior  region  in  the  equivalent  problem.  The  conductor  is  now  modeled  as  an 
impedance  sheet  at  the  conductor  surface,  as  shown  in  Fig.  lb. 


Figure  1 :  Comparison  of  the  actual  geometry  and  various  surface  equivalent  models. 

(a)  actual  geometry  with  the  descriptions  of  fields  and  material  properties;  (b)  the  equivalent  model  using  the 
effective  internal  impedance  (Eli)  where  the  exterior  fields  are  identical  to  the  original  problem,  the  conductor 
interior  is  replaced  by  the  exterior  medium,  and  an  impedance  sheet  is  defined  at  the  conductor  surface. 

The  volume  filament  technique  (Figure  2a)  has  been  widely  used  in  high  performance 
package  modeling  since  an  external  Green’s  function  can  be  used  to  calculate  mutual  inductances 
between  conductors  regardless  of  their  different  conductivities,  avoiding  solving  complicated 
boundary  condition  problems  that  are  often  encountered  in  surface  equivalent  problems.  Also, 
this  technique  directly  calculates  resistances  and  inductances  (rather  than  field  quantities)  that  are 
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frequently  required  for  circuit  simulation.  According  to  the  equivalence  theorem,  the  original 
problem  can  be  replaced  with  equivalent  surface  currents  as  long  as  this  current  generates  the 
correct  fields  outside  of  the  conductors,  while  the  fields  inside  of  the  conductor  can  be  assumed 
to  be  arbitrary.  This  allows  one  to  replace  the  medium  inside  the  conductor  boundaries  with  the 
same  medium  as  the  external  region.  Consequently,  an  external  Green’s  function  can  be  used 
everywhere  in  the  integral  equations.  ^ 

/''V'N  /H\ 


□□BDDDpED_  I  *  j 

--F — “F"F  i  . .  1 . 

Figure  2:  (a)  Illustration  of  the  volume  filament  technique  where  Ll}  represents  the  self  ( i  -  j )  and  mutual  0  *  J ) 
inductances  between  filaments;  R  is  filament  resistance  given  by  \foA  where  A  is  the  cross-sectional  area  of  the 
filament,  (b)  Surface  ribbon  method  where  L.  represents  the  self  and  mutual  inductance  between  surface  ribbons 

and  Zeii  is  the  diagonal  elements  of  [Zeii  ]  given  in  (16). 

In  formulating  a  surface  boundary  condition  the  fields  inside  of  the  conductor  are  usually 
assumed  to  be  zero  [16],  but  here  we  assume  the  vector  magnetic  potential  is  continuous  across 
conductor  surfaces,  and  the  equivalent  surface  current  is  defined  as  the  difference  between 
original  tangential  magnetic  field  and  internal  tangential  magnetic  field  generated  in  this  new 
equivalent  problem.  Figure  1  shows  the  original  problem  and  the  equivalent  problem  explained 
here.  Of  course,  it  might  seem  that  this  new  approach  would  make  problem  much  more 
complicated,  but  as  will  be  demonstrated  in  this  paper,  the  surface  impedance  boundary  condition 
satisfying  this  new  problem,  the  effective  internal  impedance  (Eli),  is  much  more  localized  than 
the  SIBC  from  low  to  high  frequencies. 

The  surface  integral  equation  can  be  obtained  from  Green’s  theorem  [3] 

JJ^’(d>V2'F-'FV2d>)=  f •  0) 

For  the  exterior  region  where  Green’s  function  satisfies  V2G=— 5{r  —  r )  and  magnetic  vector 
potential  satisfies  Laplace’s  equation,  V2  Az  =0 ,  the  equation  (1)  becomes 

||d^’[/lJ.(r’)(5(r-r’)]=-|dr’  G(r,r')  —A.(r’)  -  ,  (2) 

where  ¥  =  Az ,  0>  =  G ,  and  negative  sign  in  right  side  of  (2)  indicates  derivative  with  respect  to 
normal  pointing  inwards.  Similarly,  for  the  interior  region  where  Green’s  function  satisfies 
V 2G -  jcojucr G=-S(r- r’)  and  longitudinal  electrical  current  satisfies  V2JZ  -jco^ia  Jz  =0,  the 
equation  (1 )  becomes 


|  -  r)]= J  drfor,  • 


where  VF=  J,  =  <jE7  and  =  G . 
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The  surface  integral  equation  is  formulated  on  the  surface  of  the  conductors,  which 
requires  appropriate  distribution  of  the  delta  functions  in  the  equations  (2)  and  (3).  Suppose  that 
at  given  location  the  strength  of  the  delta  function  just  outside  of  the  surface  is  p{r').  Then  the 
strength  of  the  delta  function  just  inside  of  the  surface  is  1  -  p{r') .  The  left  side  of  the  equation 
(2)  equals  to  p{r)Az  (r),  which  can  be  expressed  as  jdr'[p(r’)  Az  {r')S(r  -  r ’)] .  After  applying 


VxA=ji  H ,  the  surface  integral  equation  for  the  exterior  region  becomes 

| dr\G\„  (r,  (r’)]-  |</r’[Ge2„  (r,  r’)-  p{r')S{r  -  r ’)]az  (r’)  =  0  ,  (4) 

where  H,  and  Az  are  the  tangential  magnetic  field  and  longitudinal  magnetic  vector  potential 
respectively,  p(r')  is  the  position  dependent  coefficient  describing  source  strength  distribution 
just  outside  of  the  conductor  surface,  and  Gl{ja  and  G]„  are  the  free  space  external  Green  s 
function  and  its  derivative  with  respect  to  the  surface  normal,  respectively,  and  the  integral  is 
carried  out  along  the  perimeter  of  all  conductors.  Similarly,  for  the  interior  region,  the  left  side 
of  the  equation  (3)  equals  to  [l-.p(r)]£z(r),  which  can  be  expressed  as 

|{[l-p(r’)]£;:(r’)<5(r-  r’)}.  After  applying  VxE=-jo)^iH  ,  the  integral  equation  for  the  inner 

r, 

region  becomes 

Uric',  (r,  AimiH,  (r')]+  U-{c,!  (r,  r')+ [ I  -  p(r’)]S(r  -  r'fe  M  =  0  •  <5) 

where  G'q  and  Gq  are  the  Green's  function  and  its  normal  derivative  for  conductor 
q  respectively,  and  the  integral  is  carried  out  along  the  each  individual  conductor  perimeter  r, . 


Equation  (4)  and  (5)  should  satisfy 

{dr’H,  =  /,  ,  (6) 

where  /  is  tb''.  total  current  carried  by  conductor  q  . 

After  ■’.pplying  the  equivalence  theorem,  the  Green’s  function  that  appears  in  equation  (5) 
can  be  replaced  by  the  free  space  Green’s  function  used  in  (4).  In  this  paper  we  use  continuity  of 
the  vector  magnetic  potential  across  conductor  surfaces  to  complete  the  definition  of  the  EH. 
Equation  (4)  for  the  exterior  region  becomes 

| dr'[G'ex, (r,  r'^H™ (r ’)]- 1 rfrfe (r,  r’) -  p{r')S(r  -  r ’)]Az  (r’)  =  0  ,  (7) 

where  Hr"  is  the  external  tangential  magnetic  field.  For  the  interior  region  where  internal 
medium  is  replaced  with  external  medium  equation  (5)  becomes 

(8) 

jdr{Gl,  (r,  r')uH"  (/)]-  §dr{G*„{r, r')+ [l  -  p{r')]5{r  -  r')K  (r')  =  0 
where  is  the  internal  tangential  magnetic  field  and  the  integral  is  now  carried  out  along  the 


entire  surface  of  inside  of  the  conductor. 

If  the  values  of  external  tangential  magnetic  field  and  magnetic  vector  potential  are 
exactly  known  in  equation  (7)  (here,  magnetic  vector  potential  is  related  to  magnetic  field  by 
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B  =  VxA),  the  position  dependent  coefficient  p{r’)  can  be  determined,  and  then  H"  can  be 
found  from  equation  (8). 

For  this  new  equivalent  problem,  the  relationship  between  the  longitudinal  electric  field 
and  the  difference  of  the  tangential  magnetic  fields  across  the  equivalent  surface  impedance  sheet 
is  called  the  effective  internal  impedance  (the  EH,  to  distinguish  it  from  the  SIBC),  and  is  then 


Mr’)= 


KM 


(9) 


where  Ztii  is  the  EH  and  Js  is  the  equivalent  surface  current  density. 

To  compute  the  integrated  quantities  of  series  impedance  using  the  EH  approach,  power 
applied,  power  dissipated,  and  magnetic  energy  stored  in  the  equivalent  problem  are  found  using 


->  ,  ,  ->  ->  noi 

=jsdsZeii\J,\  +ja)\sdsAJ's 

where  S  is  the  conductor  surface  and  O  is  the  electric  potential.  The  left-hand  side  of  the 
equation  (10)  is  the  power  applied,  the  first  term  of  right-hand  side  is  the  power  dissipated,  and 
the  second  term  corresponds  to  the  magnetic  energy  stored.  From  the  power  dissipation  term  in 
(10),  resistance  can  be  calculated  using 

l  dr'Re\Z'U\J,\2\  ,  (11) 

*r>)=- - * — 1< 

idrJs 

where  lq  is  the  length  of  the  line  q .  Internal  inductance  of  conductor  q  and  total  external 
inductance  are  calculated  from  the  stored  magnetic  energy  term. 
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where  Vexl  is  the  volume  exterior  to  the  conductors.  The  surface  inductance  due  to  the  magnetic 
energy  stored  by  the  impedance  sheet  is  given  by 

£*'i">kw!)  •  (>*> 

Z.™  (»)■=* — j— - =4— 

T«l  (O^dr'J, 

Finally,  the  sum  of  internal,  external,  and  surface  inductances  from  (12),  (13),  and  (14)  gives  the 
total  series  inductance  for  the  transmission  line  in  the  EH  formulation. 

To  validate  equations  (1 1)  to  (14),  a  similar  procedure  can  be  used  to  derive  the  equations 
for  series  impedance  using  the  boundary  element  method  [3,  17].  Using  the  exact  exterior 
magnetic  and  electric  fields  from  the  original  problem,  it  can  be  shown  that  the  resistances  given 
by  (11)  are  identical  to  the  original  problem,  that  the  external  inductance  given  by  (13)  is 
identical  to  the  original  problem,  and  the  sum  of  (12)  and  (14)  is  identical  to  the  internal 
inductances  in  the  original  problem. 

III.  Surface  Ribbon  Method 

For  series  impedance  calculation,  even  though  (11)  to  (14)  are  useful  to  validate  the  new 
approach,  it  is  a  time-consuming  task  to  solve  the  coupled  integral  equations.  A  much  simplified 
formulation  can  be  derived  directly  from  (7),  (8),  and  (9)  on  the  surface  of  the  conductor.  By 
subtracting  (8)  from  (7)  and  applying  (9)  and  £Z+V :4>=-  joAz  at  the  boundary  of  the 
conductors,  the  following  integral  equation  results  for  m  conductor  system. 


ft  dr'\jw^*S(r-ryjA}j,(rt)*'tidr'S(r-ry^'=0  .  (15) 

,=!  <7=1  * 

Each  conductor  perimeter  T,  can  be  further  divided  into  Nq  segments,  each  with  sub  section 
Cq  k ,  where  these  pieces  represent  current-carrying  “ribbons”  of  width  wqk  =  jdr .  Equation 

(15)  then  becomes 

m 

ZqJ  (r )j  *'*  (r)+  X  JA 7 ! J  (r’)Gi* (r’ r ’)  =  "V z*9*  •  (16) 

.=1  j=\  Li-' 

where  ZqJ  and  Jqk  are  respectively,  the  effective  internal  impedance  and  sheet  current  density 

on  the  k*  ribbon  (k  =  1,2, —,Nq )  of  the  q*  conductor  (q  =  1,2,..., m )  at  a  given  frequency,  and 

the  second  term  on  the  left  hand  side  represents  contributions  from  self  and  mutual  inductances. 

For  a  two  dimensional  case  G\a  (r,  r’)  =  — In|r  -  r| .  If  the  ribbons  are  narrow  enough  that  the 

L7Z 

sheet  current  density  is  constant  across  each  ribbon,  integrating  over  the  k*  ribbon  yields 
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where  /  *  is  the  total  sheet  current  carried  by  the  k*  ribbon,  and  for  simplicity  the  double 
superscript  (<7,  k )  has  been  replaced  with  the  single  superscript  k  (k  =  1,2,. . N ,  N  —  N^xm). 
Finally, 


7k  I  k  — 

" eii 1  s 


- + L  ldr'drG'<*  (r’r’)=<I>i  > 


where  Of  -Olj  is  voltage  drop  along  the  ribbon  k,  Zkn  is  the  EH  averaged  over  the  width  of 
ribbon  k,  and  lk  is  the  length  of  ribbon  k,  which  is  assumed  to  be  1  in  a  two  dimensional 
problem.  Equation  (18)  can  be  expressed  as  an  N  x  N  matrix  equation 


[fch  ;<«[/]= M  _  <19> 

where  [Zeii]  is  an  N  x  N  diagonal  matrix  made  up  of  the  Zku/wk  and  [l]  is  a  matrix  consisting 

of  self  and  mutual  inductances  between  surface  ribbons.  This  approach  is  called  the  surface 
ribbon  method  (SRM).  Now  rather  than  calculating  resistance  and  inductance  from  equation  (11) 
to  (14),  if  [Z,,,  ]  is  known  (19)  can  be  used  to  calculate  resistance  and  inductance  directly,  as  is 
done  in  the  volume  filament  method  [1,  2].  As  shown  in  the  following  section,  the  EH  has 
relatively  localized  characteristics,  allowing  [Ztii]  to  be  easily  approximated. 

Just  as  in  the  SBEM,  in  principle  it  is  as  difficult  to  find  a  priori  the  exact  EE,  especially 
at  low  and  middle  frequencies,  as  solving  the  original  problem.  However,  it  will  be  shown  below 
that  the  surface  impedance  of  a  single,  isolated  conductor  is  a  better  approximation  to  the  EH 
than  to  the  SBC  for  multi-conductor  transmission  lines  at  low  frequencies,  while  at  high 
frequencies  the  isolated  conductor  surface  impedance  is  an  adequate  approximation  of  both  the 
EH  and  the  SBC.  Hence,  more  accurate  approximation  of  the  transmission  line  series 
impedance  can  be  obtained  over  a  wider  frequency  range  using  the  more  easily  calculated 
isolated  conductor  surface  impedance  as  an  EH  in  the  SRM  than  can  be  achieved  using  lire  same 
isolated  surface  impedance  as  the  SBC  in,  for  instance,  the  SBEM. 


IV.  Sample  Results 

To  calculate  the  exact  EH,  the  volume  filament  technique  can  be  used  to  find  the  external 
magnetic  and  electric  fields  and  the  vector  magnetic  potential,  then  these  values  can  be  used  in 
(7)  and  (8)  to  calculate  the  internal  magnetic  field,  and  finally  the  EH  is  found  using  (9).  Fig.  3 
shows  a  comparison  between  the  EH  values  for  two  circular  conductors  with  radii  1  mm 
separated  by  0.2  mm,  20  mm,  and  200  mm,  and  the  surface  impedance  of  an  isolated  circular 
conductor  given  by  [18], 


where  a  is  the  radius  of  the  conductor.  The  most  interesting  characteristics  of  the  EH  is  the 
position  independence  of  the  real  part  at  low  frequency,  as  shown  in  Fig.  3a.  All  real  parts  of  the 
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FIT  are  almost  constant  and  are  quite  close  to  the  surface  impedance  of  the  isolated  conductor 
calculated  from  (20).  The  imaginary  parts  do  have  a  position  dependence,  especially  when  two 
conductors  are  vety  close  to  each  other.  However,  this  dependence  diminishes  quickly  as  the  two 
conductors  are  separated,  with  values  of  Im(EII)  being  approximately  the  average  of  that 
obtained  from  (20).  At  high  frequency,  the  internal  fields  of  the  conductors  approach  zero,  and 
the  EH  calculated  from  (9)  is  expected  to  approach  to  the  SIBC  which  is  accurately  approximated 
by  (20)  in  the  high  frequency  limit.  Fig.  3b  shows  that  this  does  indeed  occur.  In  other  words, 
(20)  can  be  used  as  an  approximation  of  the  EH  to  calculate  series  impedance  of  circular 

conductors  regardless  of  the  frequency  or  proximity  of  other  conductors  (i.e„  Z‘,  =ZCI>  in  (18) 
and  (19)).  To  support  this  point,  the  series  impedance  of  two  circular  conductors  separated  by  0.2 
mm  is  calculated  and  compared  with  the  volume  filament  method  as  shown  in  Fig.  4.  Here,  a 
polygon  with  degree  24  shown  in  Fig.  4a  is  used  instead  of  a  circular  conductor  to  simplify  the 
calculation  of  the  mutual  inductance  between  ribbons.  Figure  4b  indicates  (20)  is  indeed  a  good 
approximation  of  the  EH  for  series  impedance  calculation  even  when  conductors  are  very  close  to 
each  other. 


Figure  3(a)  Figure  3(b) 

Figure  3:  Position  dependency  of  Eli  for  circular  conductor  (copper)  with  a  =  1mm,  d  =  0.2, 20,  and  200  mm.  solid 
line:  d  =  0.2  mm,  dotted  line:  d  =  2  mm,  dashed  line:  d  =  200  mm,  dot-dash  line:  equation  (17).  (a):  low  frequency 
limit  (1  kHz);  (b):  high  frequency  limit  (1  MHz). 
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Figure  4(a)  '  Figure  4(b) 

Figure  4:  (a)  Segmentation  scheme  for  circular  conductor  using  polygon  of  degree  24  in  the  SRM.  (b)  Series 
impedance  calculation  of  two  closely  spaced  circular  conductors  (copper)  with  a  =  1  mm,  d  =  0.2  mm,  N  =  48.  Solid 
line:  SRM;  dashed  line:  VFM 

The  EH  values  for  square  conductors  with  width  25  pm  and  separated  by  5  pm  and  50  pm 
have  also  been  calculated.  It  is  impossible  to  calculate  exactly  the  SIBC  of  isolated  square  or 
rectangular  conductors  due  to  the  singularity  at  the  comer,  however,  it  can  be  approximated  [19- 
21].  Fig.  5  shows  the  low  frequency  characteristics  of  the  EH.  Both  real  and  imaginary  parts  of 
the  EH  show  similar  characteristics  as  the  circular  case  except  near  the  comer. 
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Figure  5(a)  Figure  5(b) 

Figure  5:  Position  dependency  of  Eli  for  square  conductor  (copper)  with  25  pm  width  and  d  =  5  pm  and  50  pm  at  f = 
1  kHz.  solid  line:  d  =  5  pm,  dashed  line:  d  =  50  pm.  (a)  real  part  (b)  imaginary  part 

At  high  frequency,  the  surface  impedance  of  an  infinitely  wide  and  thick  conductor  given 
by 


is  a  reasonable  approximation  for  the  SIBC  far  from  comers;  Fig.  6  shows  that  the  EH  also 
approaches  (21)  except  near  the  comer.  Examination  of  other  geometrical  cases  has  confirmed 
that  the  EH  is  well  approximated  by  the  surface  impedance  of  an  isolated  conductor,  regardless  of 
the  presence  of  other  conductors.  Various  studies  [17, 22]  have  shown  that  an  approximated  EH 
can  be  directly  applied  to  (19)  to  accurately  estimate  the  series  impedance  of  various  cross- 
sectional  conductors.  The  series  impedance  calculations  for  multi-conductor  transmission  lines 
using  the  SRM  are  detailed  in  [17, 22]. 
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Position  s  (meters)  Position  s  (meters) 

(a)  (b) 

Figure  6:  Position  dependency  of  Eli  for  square  conductor  (copper)  with  25  pm  width  and  d  =  5  pm  and  50  pm  at  f  - 
1  GHz.  solid  line:  d  =  5  pm,  dotted  line:  d  =  50  pm,  dashed  line:  equation  (18).  (a)  real  part  (b)  imaginary  part. 


V.  Conclusions 

The  precise  definition  of  an  alternative  surface  impedance  boundary  condition  has  been 
formulated  for  use  in  lossy  transmission  line  analysis.  For  this  boundary  condition,  the  conductor 
is  modeled  as  an  impedance  sheet  on  the  conductor  perimeter,  the  conductor  interior  is  replaced 
by  a  (non-conducting)  exterior  material,  and  the  effective  internal  impedance  (EH)  is  defined  on 
the  impedance  sheet.  At  low  frequency,  the  real  part  of  the  EH  is  almost  constant  along  the 
perimeter  of  a  conductor,  and  at  high  frequency  the  EH  approaches  the  SIBC,  hence,  the  EH  is 
approximately  a  local  quantity  from  low  to  high  frequency.  Approximations  to  the  EH  can  be 
used  in  the  surface  ribbon  method,  producing  numerically  efficient  and  accurate  calculation  of 
the  series  impedance  of  lossy  multi-conductor  lines.  Here,  die  SRM  was  applied  to  two 
dimensional  problems,  but  can  be  extended  to  three  dimensional  problems  [17].  Since  this 
method  requires  substantially  less  memory  than  the  VFM,  the  SRM  can  be  used  effectively  for 
ground  plane  modeling  including  current  crowding  effects. 

Table  1  compares  the  number  of  unknowns  and  run  time  on  an  IBM  RISC  6000  for  the 
FBEM,  the  I-BEM,  and  the  I-SRM.  The  surface  techniques  are  at  least  25  times  faster  than  the 
FBEM  in  assembling  the  required  matrices,  and  at  least  11  times  faster  in  solving  using  a 
gaussian  elimination  algorithm.  In  addition,  the  SRM  and  the  SBEM  require  matrix  assembly 
only  once  for  all  frequencies  of  interest,  while  FBEM  requires  assembling  the  matrix  at  each 
frequency.  Hence,  the  surface  techniques  are  significantly  more  computationally  efficient,  and 
by  using  the  isolated  conductor  surface  impedance  as  an  approximation  for  the  EH  in  the  SRM, 
little  accuracy  is  lost  over  the  whole  low-to-high  frequency  range. 
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Method 

Number  of 
unknowns 

CPU  time  (sec) 
Assembling  Solving 

BEM 

578 

468.1 

100.10 

I-SBEM 

290 

18.3 

I-SRM 

287 

5.1 

8.96 

Table  1 :  Run  time  comparison  between  BHM  and  SRM. 

Comparison  of  run  time  on  an  IBM  RISC  6000  for  BEM,  BEM  assuming  the  SIBC  is  equal  to  the  isolated  conductor 
surface  impedance  (1-SBEM),  and  SRM  assuming  the  Eli  is  equal  to  the  isolated  conductor  surface  impedance  (I- 
SRM).  The  surface  of  a  conductor  is  segmented  into  144  elements.  I-SBEM  and  I-SRM  are  at  least  25  times  faster 
than  BEM  in  assembling  a  matrix,  and  at  least  11  times  in  solving  the  matrix  using  a  simple  gaussian  elimination 
algorithm. 

In  summary,  this  work  has  shown  a  new  surface  impedance  technique  that  can  be 
accurately  approximated  using  the  isolated  conductor  surface  impedance,  even  at  low  frequency. 
This  effective  internal  impedance  (EH)  can  be  combined  with  the  surface  ribbon  method  (SRM) 
allowing  numerically  efficient  and  accurate  calculation  of  the  series  impedance  of  lossy  multi¬ 
conductor  transmission  lines.  In  this  technique,  the  effective  internal  impedance  defines  an 
impedance  sheet  at  the  conductor  surfaces  and  the  conductor  interiors  are  replaced  with  the 
exterior  dielectric.  For  circular  cross  section  conductors,  using  the  surface  impedance  of  an 
isolated  circular  conductor  in  the  EC-based  SRM  yields  an  excellent  approximation  from  low  to 
high  frequency.  This  technique  can  also  be  easily  applied  to  any  lossy  multi-conductor  line 
structure  using  polygonal  cross-section  conductors  by  using  an  appropriate  effective  internal 
impedance  model. 
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SIMIAN:  a  Two  Dimensional  Multi-Conductor  Interconnect  Series  Impedance 
Calculator 

SIMIAN  (Surface  Impedance  Method  for  Interconnect  Analysis)  is  a  two  dimensional 
frequency  dependent  series  impedance  extraction  tool  for  interconnects  and  transmission  lines 
using  conductors  of  rectangular  cross  section  based  on  the  EH  and  SRM  discussed  above.  Unlike 
the  Volume  Filament  Method  [19],  only  the  surface  of  the  conductor  is  divided,  which  can  save 
significant  amounts  of  computation  time,  especially  when  the  cross  sectional  size  of  the 
conductor(s)  is  (are)  comparable  to  the  skin  depth  [14].  This  method  has  also  been  extended  to 
solve  three  dimensional  problems  [15].  SIMIAN  returns  the  frequency  dependent  series 
impedance  matrix  for  an  n-conductor  system,  i.e.,  the  n  x  n  matrix  consisting  of  self  and  mutual 
resistances  and  self  and  mutual  inductances  of  and  between  all  n  conductors. 

The  SIMIAN  program  was  successfully  complied  with  ANSI  C(g++).  The  platforms 
tested  thus  far  include: 

•  Sparc  20  machine  with  Sunos  4.1.3 

•  Sparc  20  machine  with  Solaris  2.5 

•  Pentium  machine  with  Linux 

•  HP  workstation 

The  progran  can  be  obtained  over  the  World  Wide  Web  at  the  url: 

http://weewave.mer.utexas.edu/MED_files/MED_research/Intrcncts/SIMIAN_stuff/simian_links. 

html 

Since  the  total  number  of  ribbons  used  for  each  conductor  determines  the  computation 
time,  it  is  very  important  to  reduce  them  without  loss  of  accuracy.  Since  SIMIAN  only  divides 
the  conductor  on  its  surface,  the  problem  can  be  reduced  from  N2  to  4N  compared  to  the  volume 
filament  method.  However,  more  significant  advantage  comes  from  reduction  of  N  itself.  To  get 
accurate  result  in  the  volume  filament  method,  each  segment  has  to  be  comparable  in  size  to  the 
skin  depth  at  the  frequency  of  interest;  if  a  filament  is  much  larger  than  about  1/3  of  a  skin  depth, 
serious  error  can  result.  However,  in  the  surface  ribbon  method  a  much  smaller  number  of 
segments  can  be  used  while  maintaining  reasonable  accuracy  (typically  less  than  1%  error  at 
“low”  frequencies,  perhaps  as  much  as  10%  over  a  narrow  band  of  frequencies  at  “mid” 
frequency,  and  again  about  1%  at  “high”  frequency).  For  most  cases  using  more  than  five 
ribbons  on  each  side  would  not  be  necessary,  and  in  many  cases  as  few  as  one  per  face  is 
necessary. 

Figure  7  shows  minimum  segmentation  method  [20]  using  ‘plate’  option  in  Conductor 
type  variable  applied  to  ground  plane.  Only  one  ribbon  is  used  on  each  side  of  signal  line 
resulting  in  a  total  of  20  ribbons. 


Figure  7:  Minimum  Segmentation  Scheme  with  ground  plane. 
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Examples  and  Results 
Square  Twin  Lead 

Example  1  is  simulated  with  SIMIAN  and  compared  with  FastHeniy  [21]  that  uses  the  volume 
filament  method  accelerated  using  a  multi-pole  algorithm.  As  Fig.  8  shows,  SIMIAN  can  reduce 
matrix  si?**  significantly.  For  the  problem  given,  SIMIAN  only  uses  30  to  40  total  ribbons  while 
the  volume  filament  method  needs  more  than  800  filaments  to  accurately  predict  total  resistance 
at  high  frequency.  FastHenry  accelerates  matrix  calculation  with  the  help  of  a  multi-pole 
algorithm,  and  this  same  algorithm  can  also  be  applied  to  the  surface  ribbon  method.  Figure  9 
shows  the  results  of  FastHenry,  SIMIAN,  and  the  surface  ribbon  method  (SRM)  using  a  multi¬ 
pole  algorithm.  A  total  40  ribbons  were  used  for  SIMIAN  and  SRM  with  multi-pole  algorithm, 
while  800  filaments  were  used  for  FastHenry.  In  this  particular  example,  the  time  advantage  can 
be  400  to  8000  times  faster  than  using  the  volume  filament  method,  depending  on  what  method 
is  used  for  matrix  calculation.  Considering  this  is  for  only  one  frequency  point,  tremendous  time 
savings  can  be  achieved  by  using  the  SRM. 


Figure  8:  (a)  FastHenry  [21]  results  versus  number  of  filaments;  (b)  SIMIAN  versus  number  of  ribbons. 
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Figure  9:  solid  line:  SIMIAN  (40  ribbons);  circle:  SRM  with  multi-pole  algorithm  (40  ribbons);  cross:  FastHenry 
(800  filaments),  (a)  Real  part  of  total  Impedance  (R).  (b)  Imaginary  part  of  total  Impedance  (tflL). 


Three  Conductors  over  Finite  Size/Conductivity  Ground  Plane 

Figure  10  shows  results  from  simulating  the  three  conductor  problem  shown  below.  The  total 
number  of  ribbons  used  in  SIMIAN  was  20,  which  was  considerably  less  than  FastHenry  (1250). 
Reasonable  accuracy  was  achieved  through  SIMIAN  within  a  second  (20  frequency  points  with 
SPARC  20). 
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R12(0hm/m)  R11  (Ohm/m) 


Hz 


(a)Example  2,  Z1 1 
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(c)Example  2,  Z13 


Figure  10:  Real  and  Imaginary  part  of  Series  Impedance  of  Example  2. 
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Efficient  Time  Domain  Boundary  Conditions  for  Finite  Conductivity  Metals 

As  clock  speeds  become  faster  and  chips  become  more  complex,  it  is  more  difficult  to 
accurately  predict  delay  and  crosstalk  in  densely  packed  multi-conductor  interconnects.  The 
impact  of  frequency  dependent  inductance  and  resistance  due  to  finite  conductivity  conductors  on 
time  domain  waveforms  has  often  been  ignored  or  only  poorly  approximated  even  though  this 
can  cause  overestimation  of  far-end  crosstalk  and  underestimation  of  delay.  To  include  frequency 
dependent  effects  accurately  in  time  domain  calculations,  transformation  using  FFT  is  usually 
required.  However,  it  is  difficult  to  implement  non-linear  drivers  using  this  method.  In  [3],  a 
surface  impedance  boundary  condition  was  applied  directly  to  the  time  domain  to  solve  radiation 
problems  where  high  frequency  approximation  for  surface  impedance  boundary  condition  is 
accurate  enough.  Here  we  show  a  new  method  for  transmission  line  analysis  where  the  boundary 
condition  is  represented  directly  in  the  time  domain.  Unlike  other  time  domain  simulation 
methods  in  which  all  the  RLC  parameters  of  the  transmission  lines  must  be  provided  a  priori,  in 
our  model  all  required  values  are  determined  directly  and  simply  from  conductor  geometry.  The 
resulting  frequency  domain  model  is  valid  over  a  wide  frequency  range,  and  hence  giving 
accurate  time  domain  results.  The  method  is  also  compatible  with  non-linear  drivers  and  loads. 

Application  of  impedance  boundary  condition  in  the  time  domain 

Rapid  calculation  of  interconnect  series  impedance  in  the  frequency  domain  can  be 
performed  using  an  impedance  boundary  condition  [2].  For  a  rectangular  conductor  cross  section, 
the  cross  section  is  partitioned  into  two  trapezoidal  and  triangular  parts  as  shown  in  figure  2a;  the 
“transverse  impedance”  of  each  part  is  used  as  an  approximation  of  impedance  boundary 
condition  for  the  conductor  [2].  To  complete  the  problem,  each  transverse  impedance  is  assigned 
to  a  corresponding  surface  ribbon  and  mutual  inductances  between  these  ribbons  are  calculated. 
The  resulting  two  dimensional  equation  is 

•  <■> 

where  [Zeii]  is  a  diagonal  matrix  of  the  effective  internal  impedances  (EH)  (representing  the 
impedance  boundary  condition)  and  [L]  is  a  dense  matrix  of  mutual  inductances  between 
ribbons. 

To  use  this  frequency  domain  concept  in  time  domain  analysis,  the  equivalent  circuit 
shown  in  figure  2b  is  used  to  model  each  ribbon.  This  produces  a  rational  function  in  the  s- 
domain  that  is  easily  converted  into  time  domain  exponential  functions.  All  values  in  this  circuit 
model  are  determined  directly  from  simple,  closed  form  expressions  depending  on  only  the 
conductivity  and  geometry  of  the  conductors  [5, 6],  Transforming  (1)  into  the  time  domain  gives 

where  **’  represents  convolution,  [^(f)]  is  a  sum  of  exponential  functions  obtained  by  converting 
the  rational  function  in  s-domain  [Ze///s],  i.e.. 
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(3) 


Z,.y  _  a3*3+<V2+<M  +  flo  -.*■  |  *2  ,  Jh—+Ji i—  . 

J  i(V3  +b2S2  +fc,j  +  fe0)  5  •S  +  /72  *  +  F3  ^+P4 

The  rational  function  in  (3)  can  be  further  reduced  by  using  dominant  pole  approximation  to 
reduce  the  computation  time  and  memory  usage. 

The  convolution  appearing  in  (2)  can  be  evaluated  easily  using  a  recursive  relationship 

m, 

Y(nAt)  =  kepAI  ■  X(nAt )  =  kAt- X(nAt)+ep&tY((n-1)At)  .  (4) 

Applying  (3)  to  (2)  yields 

(JM(]+[L])|  [/]+[^s]=  [i]|  [/MM=  M  ,  (5) 

where  (Vds  ]  is  an  extra  voltage-dependent  source  added  to  the  telegrapher’s  equation,  depending 

on  Af  and  on  the  poles  and  residues  of  the  equivalent  circuit  of  the  EH.  Since  At  is  usually  very 
small,  this  voltage  dependent  source  is  mainly  responsible  for  incorporating  the  current 
redistribution  effects  (i.e.,  the  frequency  dependent  proximity  and  skin  effects)  in  the  time 
domain.  Equation  4  is  then  solved  along  with  the  second  telegrapher’s  equation  for  the  shunt 
admittance  (e.g.,  the  capacitance)  of  the  interconnect. 

Method  of  characteristics 

The  method  of  characteristics  was  first  developed  to  simulate  a  single  lossless  line  [8], 
then  later  extended  to  multi-resistive  lines  [9].  Recently,  it  was  further  extended  to  non-uniform 
lossy  line  cases  [10],  where  frequency  dependent  series  impedance  is  calculated  using  FFT. 
Here,  for  simplicity,  a  homogeneous  medium  with  [G]=0  is  assumed.  More  general  cases  are 
well  demonstrated  in  [10].  From  equation  4  and  the  corresponding  admittance  equation, 
d([V]±[Rll]ydt  =  -[Vds}vp  ,  (6) 

where  [R]=vp\L]vp  is  propagation  velocity  and  +/-  represent  incident  and  reflected  waves, 

respectively.  To  capture  the  distributed  resistive  effects  intermediate  segmentation  is  required, 
then  equation  8  becomes 

for  the  incident  wave,  and  for  the  reflected  wave 

•IV]  -  =  {[V  ]-  PI']  +  Af  l/p  [Kk  ] f„  )  .  (7b) 

The  equivalent  model  is  shown  in  figure  3;  comparing  to  the  model  shown  in  [11],  the  resistance 
is  replaced  with  a  new  voltage-dependent  source  to  capture  frequency  dependent  effects.  To 
simulate  rectangular  conductors,  each  conductor  is  divided  into  four  parts  as  shown  in  figure  2a, 
then  each  part  is  represented  with  its  equivalent  circuit,  and  finally  four  parts  are  connected  in 
parallel. 
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Results  and  conclusions 

Figure  4  shows  simulation  results  obtained  using  the  new  method  introduced  in  this 
paper,  as  well  as  illustrated  the  errors  resulting  from  failure  to  properly  include  the 
time/frequency  dependencies  of  both  resistance  and  inductance.  Excellent  agreement  is  shown 
compared  to  FFT  method  using  a  full  dispersion  curve,  while  the  CPU  time  required  for  this 
method  (7.6  seconds)  is  much  less  than  the  FFT  method  (283  seconds).  In  conclusion,  this  paper 
demonstrates  that  the  dc  to  skin  effect  frequency  variation  of  series  impedance  can  be  directly 
and  accurately  included  in  the  time  domain,  with  all  needed  values  can  be  easily  determined 
directly  from  conductor  geometry,  while  maintaining  high  numerical  efficiency. 


Figure  2)  (a)  Segmentation  scheme  for  rectangular  bar  to  approximate  surface  impedance,  (b) 
Equivalent  circuit  model  for  approximated  surface  impedance. 


R  R  Vk  r  R 


Figure  3)  Equivalent  model  for  method  of  characteristics.  ER  and  £/  are  voltage  dependent 
sources  of  reflected  and  incident  waves  respectively  and  appear  in  second  part  of  right  hand  side 
of  the  equation  8a  and  8b,  respectively. 


AFOSR  Grant  No.  F49620-96- 1-0032 


23 


Final  Report 


(a)  far  end  waveform  at  line  1  and  2;  solid 
line:  FFT  result,  dotted  line:  new  method 
described  here. 
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line:  FFT  result,  dotted  line:  new  method 
described  here. 


(c)  far  end  voltage  waveform  at  line  4:  errors 
induced  by  failure  to  correctly  include  both 
resistance  and  inductance  frequency 
dependencies:  Solid  line:  FFT  using  full 
dispersion  curve;  dotted  line:  constant  R-L-C 
transmission  line  analysis;  dashed-dotted 
line:  Star-Hspice®  97.2  using  W  Element 
where  one  high  frequency  resistance  value 
must  be  calculated  using  some  other  method 
(here,  the  surface  ribbon  method  was  used 
[2]). 


Figure  4)  Simulation  results  for  4  lines  above  a  ground  plane.  Conductor  geometry:  width  &  thickness 
=  20  pm,  separation  between  conductors  =  1 0  (im,  signal  conductors  10  pm  above  ground  plane, 
ground  plane  20  pm  thick,  170  pm  wide,  line  length  =  0.1  m.  Input  wave  form  applied  to  line  1  only: 
IV  trapezoidal  wave  with  rise,  fall  time  =  0.1  ns,  duration  =  1  ns;  source  resistance  =  5  0,  terminated 
with  10  pF  capacitance. 
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HP-Finite  Elements  Applied  to  Electromagnetic  Problems 

Another  important  advance  in  the  area  of  computational  electromagnetics  for  application  to 
mixed  signal  simulation  has  been  in  the  area  the  finite  element  method  (FEM).  A  new  hp-adaptive 
technique  Ivas  recently  been  validated  that  allows  the  local  variation  of  both  the  element  size  (h)  and  the 
order  of  approximation  (p).  This  approach  has  the  advantage  that  complex  structures  should  be  easily 
simulated.  The  formulation  has  been  demonstrated  to  allow  discontinuous  changes  in  material 
properties  (such  as  the  change  from  a  dielectric  to  a  metal),  to  be  stable,  and  to  easily  handle 
curvilinear  geometries.  Extensive  numerical  testing  is  still  required  in  this  area. 

Numerical  schemes  based  on  nodal  finite  elements  have  been  successfully  applied  to  many 
electromagnetic  problems  formulated  in  terms  of  potentials.  Unfortunately,  this  methodology 
frequently  fails  if  extended  directly  to  the  field  equations,  and  new  families  of  so  called  “edge-based” 
elements  have  been  proposed.  However,  in  the  low  frequency  range,  the  traditional  formulation 
becomes  unstable  even  with  the  “edge-based”  elements  and  can  lead  to  spurious  solutions.  In  this 
thesis,  we  present  a  general  hp-adaptive  methodology,  free  of  spurious  solutions  over  the  entire  range 
of  physically  relevant  frequencies  for  steady-state  Maxwell’s  equations  in  interior  domains. 

We  design  our  family  of  H(curi)-conforming  hp-edge  elements  in  the  context  of  a  model  problem  with 
sources.  These  elements  generalize  the  well-known  edge  elements  of  Nedelec  [4,  5].  We  confirm  that 
finite  element  spaces  built  on  our  elements  ensure  discrete  compactness  and  commutativity  of  de  Rham 
complex,  two  fundamental  properties  necessary  for  stability  and  convergence  of  mixed  methods  for 
Maxwell’s  equations.  Analytical  results  suggest  that  the  method  remains  stable  even  at  very  low 
frequencies. 

Following  the  logic  applied  to  the  model  problem,  we  develop  a  new  formulation  for  the  full-wave 
analysis  of  waveguides.  Using  the  discrete  compactness,  we  prove  that  this  formulation  leads  to 
numerical  approximations  which  define  a  family  of  collectively  compact  operators.  Relying  on  the 
result  by  Chatelin  [1],  we  show  that  these  operators  converge  not  only  pointwise  but  also  in  the 
operator  norm.  To  the  best  of  our  knowledge,  this  result  as  well  as  a  study  on  spectral  convergence  at  a 
given  frequency,  has  never  been  published  in  the  literature. 

In  our  study  on  spectral  convergence,  we  relate  the  discrepancy  between  exact  and  numerical 
eigenvalues  to  the  interpolation  error.  The  actual  convergence  rate  estimates  are  for  shape-regular 
meshes  with  element-dependent  order  of  polynomial  approximation  p  which  remains  uniformly 
bounded  from  above.  The  convergence  estimates  are  derived  for  non-conducting  Lipschitz  polygons 
with  piece-wise  continuous  e,  m,  and  for  polygons  with  slits  provided  e,  m,  are  constants.  We  stress 
that  our  approach  is  general  and  applies  to  any  family  of  vector  elements  which  satisfy  the  discrete 
compactness,  approximability  and  some  compatibility  properties.  For  example,  our  analysis  covers 
standard  Nedelec  elements  of  [4,  5].  Based  on  the  new  formulation,  a  finite  element  code  for  the  full- 
wave  analysis  of  waveguides  has  been  written.  This  Fortran  90  code  is  built  on  top  of  the  2Dhp90_EM 
package  for  electromagnetics  [6]  which  supports  the  new  hp-edge  finite  elements. 


AFOSR  Grant  No.  F49620-96- 1-0032 


26 


Final  Report 


Performed  numerical  experiments  fully  confirm  our  analytical  convergence  estimates  obtained  for 
affine  non-conducting  domains.  Figure  5  and  Figure  6,  and  suggest  that  these  estimates  remain  valid 
for  domains  with  nonzero  conductivity  and  curved  boundaries.  Figure  7  and  Figure  8. 

We  show  that  the  code  correctly  captures  and  tracks  the  development  of  a  complex  mode  in  a  shielded 
lossless  microstrip  in  Figure  1.  The  dispersion  curves  in  Figure  2  for  some  representative  modes, 
including  a  complex  mode,  are  produced  using  non-uniform  hp-meshes  with  two  levels  of 
h-refinement  performed  at  the  tip.  Scalar  hp-quads  of  order  ps  =  2,  3, 4  are  linked  with  compatible  hp- 
edge  elements. 

To  test  the  method  on  structures  of  practical  interest,  we  perform  the  full  wave  analysis  of  microstrip 
depicted  in  Figure  3.  Dimensionless  numerical  frequency  corresponding  to  the  domain  used  to  compute 
attenuation  constants  of  this  microstrip  can  drop  below  10-6  on  the  global  level,  and  below  10-8  on  the 
element  level.  The  obtained  results  show  that  the  formulation  is  indeed  stable  as  w  — F  0.  In  Figure  4, 
we  observe  that  the  FE-method  based  on  our  formulation  and  implemented  using  the  hp-elements 
yields  the  attenuation  curve  marked  “FE”  which  is  in  a  good  agreement  with  the  curve  computed  by 
means  of  a  conformal  mapping  technique  taken  from  [3]. 

We  believe  that  our  work  has  prepared  a  solid  foundation  for  a  fully  automatic  hp-adaptive  finite 
element  code  for  full-wave  analysis  of  waveguides.  For  a  detailed  discussion  on  the  subject,  we  refer  to 
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Experimental  Measurements  on  a  Prototypical  Personal  Communications  Services  (PCS) 
System 

Our  research  activities  in  this  area  have  been  focused  on  building  RF  transceivers  and  baseband 
circuits  for  a  handset  testbed  at  1.8GHz,  the  personal  communication  services(PCS)  band.  One 
important  objective  of  this  effort  is  to  fully  understand  how  the  PCB  layout,  chip  placement  and 
packaging  affect  the  overall  wireless  system  specifications.  This  knowledge  will  help  us  to  properly 
design  multi-chip  module  (MCM)  package  for  wireless  communication  systems. 

The  stability  of  frequency  synthesizers,  realized  by  phase  locked  loop,  is  critically  important  to 
the  quality  of  wireless  communications.  If  they  lose  lock  even  for  a  short  period  of  time,  they  can  cause 
significant  problems  with  the  modulation  and  demodulation.  In  our  experiments,  we  found  that  the 
phase  locked  loop  (PLL)  was  very  sensitive  to  the  interference  from  the  power  supply,  sudden  changes 
of  impedance,  and  interference  from  power  amplifiers.  The  RF  transceivers  we  built  were  for  the  time- 
division-duplex  (TDD)  applications,  in  which  the  transceivers  transmit  in  a  5ms  slot  and  receive  in 
another  5ms  slot  During  the  transition  period,  the  transmitter  mixer  is  turned  off  and  on  to  minimize 
the  RF  coupling  between  the  transmit  and  receive  parts.  However,  this  causes  sudden  changes  of  the 
load  impedance  to  the  PLL  and  unlocks  the  PLL.  Although  the  PLL  can  recover  from  this  sudden 
change  eventually,  it  may  cause  significant  distortion  or  phase  noise  to  the  beginning  part  of  the  signal. 
Hence,  in  a  mixed  signal  package,  it  is  better  to  build  some  driver  on  the  die  or  inside  the  package  to 
avoid  the  detrimental  effect  due  to  the  change  of  impedance.  Also,  the  PLL  must  not  share  the  same 
power  supply  as  the  power  amplifier  since  the  supply  current  surges  after  the  power  amplifier  is  turned 
on  and  dips  if  the  power  amplifier  is  turned  off.  This  is  another  source  for  the  PLL  to  lose  lock  for  a 
short  period  of  time.  Hence,  in  the  MCM  package,  we  should  have  a  different  VCC  pin  just  for  the  PLL 
power  supply.  Another  source  of  the  PLL  problem  is  due  to  the  leakage  power  from  the  power 
amplifier.  If  we  lay  the  RF  board  on  a  table,  there  is  no  problem.  However,  once  we  place  another 
baseband  board  on  the  top  of  that,  the  leaked  power  will  be  reflected  from  the  baseband  board  and 
coupled  to  the  PLL.  If  we  wish  to  package  a  PLL  in  the  same  MCM  package  as  the  power  amplifier, 
care  must  be  taken  to  isolate  the  PLL  from  the  power  amplifier. 

On  the  baseband  board,  we  used  the  Motorola  digital  signal  processor  MC56167  as  a  vocoder 
DSP.  There  are  several  delicate  design  issues  and  test  issues  we  would  like  to  address.  This  part  is  a 
mixed  signal  part.  Although  the  analog  ground  and  digital  ground  are  separate  in  the  chipset,  care  must 
be  taken  to  eliminate  the  digital  noise  in  the  analog  codec  part.  We  have  tried  numerous  PCB  layouts 
and  still  could  not  eliminate  the  digital  background  noise.  Later  we  managed  to  minimize  the  noise  by 
balance  the  load  of  the  DSP.  In  our  vocoder,  the  encoding  part  takes  much  more  time  than  the  decoder 
part.  Thus,  the  DSP  may  be  busy  for  some  time  and  idle  for  some  time.  In  this  case,  the  noise  from  the 
DSP  shows  some  periodical  pattern  and  is  more  easily  sensed  by  a  human  ear.  We  managed  to  change 
the  DSP  structure  and  made  it  busy  all  the  time.  This  makes  the  noise  white  and  not  so  sensible. 
Therefore,  for  testing  the  digital  interference  to  the  analog  circuits  in  a  mixed  signal  package,  we 
should  try  to  have  several  test  patterns  with  different  loading  to  assess  the  digital  interference  to  the 
analog  part. 
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In  summary,  we  have  extended  our  ability  to  efficiently  simulate  the  behavior  of  realistic 
interconnects  that  make  use  of  finite  conductivity  metals.  In  particular,  we  have  developed  and  tested: 

•  a  fundamental  electromagnetic  formulation  of  the  effective  internal  impedance  boundary 
condition  for  finite  conductivity  metals. 

•  an  extremely  efficient  multi-conductor  simulation  program  using  the  surface  ribbon  method  for 
interconnect  series  impedance  calculations. 

•  initial  computational  electromagnetics  code  for  FEM  simulations. 

•  experimental  work  on  mixed  signal  PCS  prototypes  to  investigate  analog  system  corruption  due 
to  cross  talk  and  power  supply  noise. 
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